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High-energy neutrino astronomy will probe the working of the most violent phenomena in the 
Universe. The Giant Radio Array for Neutrino Detection (GRAND) project consists of an array 
of ~ 10 5 radio antennas deployed over ~ 200 000 km 7 in a mountainous site. It aims at detecting 
high-energy neutrinos via the measurement of air showers induced by the decay in the atmosphere 
of T leptons produced by the interaction of the cosmic neutrinos under the Earth surface. Our 
objective with GRAND is to reach a neutrino sensitivity of 3 x lCD 11 /? -2 GeV 1 cm -2 s 1 sr 1 
above 3 x 10 16 eV. This sensitivity ensures the detection of cosmogenic neutrinos in the most 
pessimistic source models, and about 100 events per year are expected for the standard models. 
GRAND would also probe the neutrino signals produced at the potential sources of UHECRs. 

We show how our preliminary design should enable us to reach our sensitivity goals, and present 
the experimental characteristics. We assess the possibility to adapt GRAND to other astrophysical 
radio measurements. We discuss in this token the technological options for the detector and the 
steps to be taken to achieve the GRAND project. 
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Neutrinos are unique messengers that let us see deeper in objects, further in distance, and 
pinpoint the exact location of their production. These characteristics will help unveil the mys¬ 
tery behind the most violent astrophysical phenomena from which high-energy neutrinos should 
originate, and identify the source of their parent ultrahigh energy cosmic rays (UHECRs). 

The Giant Radio Array for Neutrino Detection (GRAND) aims primarily at reaching a sen¬ 
sitivity that ensures the detection of such neutrinos, even in the most pessimistic UHECR source 
models. For standard source models, it should detect enough events to launch neutrino astronomy. 

The project consists of an array of ~ 10 5 radio antennas deployed over an area of ~ 200 000 
km 2 in a mountainous site. GRAND will search for the radio signal emitted by the air showers of 
x leptons produced by the interaction of cosmic neutrinos underground. We present the detection 
method in sec. 1, the simulations that have led to the current design of GRAND in sec. 2, and 
discuss in secs. 3-4 the steps to be taken to fully define our science case, elaborated in secs. 5-6. 

1. Detection Method 

Cosmic v T s can produce x particles underground through charged current interaction, ts travel 
to the surface of the Earth and decay in the atmosphere, generating Earth-skimming extensive air 
showers (EAS) 1 [1,2]. Coherent electromagnetic radiation is associated to the shower development 
at frequencies of a few to hundreds of MHz at a detectable level for showers with E > 10 17 eV. 

The strong beaming of the electromagnetic emission combined with the transparency of the 
atmosphere to radio waves will allow the radio-detection of EAS initiated by x decays at distances 
up to several tens of kilometers. Radio antennas are thus ideal instruments for this purpose. Fur¬ 
thermore, they offer practical advantages (limited unit cost, easiness of deployment, ...) that allow 
the deployment of an array over very large areas, as required by the expected low neutrino rate. 

Remote sites, with low electromagnetic background, should obviously be considered for the 
array location. In addition, mountain ranges are preferred, first because they offer an additional 
target for the neutrinos, and also because mountain slopes are better suited to the detection of 
horizontal showers compared to flat areas, parallel to the showers trajectories. 

GRAND antennas are foreseen to operate in the 30 — 100 MHz frequency band. Short wave 
background prevents detection below this range, and above the coherence of the geomagnetic emis¬ 
sion fades. However, an extension of the antenna response up to 200 — 300 MHz would enable us to 
observe the Cherenkov ring associated with the air shower, that represents a sizable fraction of the 
total electromagnetic signal and may provide an unambiguous signature for background rejection. 

2. GRAND layout and neutrino sensitivity 

We present here a preliminary evaluation of the potential of GRAND for the detection of 
cosmic neutrinos. We initially simulated the response of a setup of 40 000 antennas deployed on a 
square layout of 60000 km 2 in a remote mountainous area (Tianshan mountains, Xinjiang, China). 

Simulation method. We perform a ID tracking of a primary v T , simulated down to the con¬ 
verted x decay. Standard rock with a density of 2.65 g/cnr 3 is assumed down to sea level and the 

'Other neutrino flavors can be neglected as the electron range in matter at these energies is too short and the muon 
decay length too large compared to flight distances underground. 
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Earth core is modeled with the Preliminary Reference Earth Model [3]. The simulation of the Deep 
Inelastic Scattering of the neutrinos is performed with Pythia6.4, using the CTEQ5d probability 
distribution functions (PDF) combined with [4] for cross section calculations. The propagation of 
the produced t is simulated using randomized values from parameterisations of GEANT4.9 PDFs 
for T path length and proper time. Photonuclear interactions in GEANT4.9 have been extended 
above PeV energies following [5]. The T decay is simulated using the TAUOLA package. 

Recently EAS radio signals have been measured in good agreement with the predictions from 
several independent simulation codes [6], but additional work is still required to properly model 
horizontal and upward-going trajectories. For the time being only an analytical parametrization of 
EAS radio detection has been performed, based on experimental results and simulation outputs. 
First, we assume that the EAS electromagnetic signal strength £ does not vary significantly along 
the shower axis for distances between 20 and 120 km from the t decay point. This is supported by 
a simulation earned out with the SE1FAS code [7] for horizontal showers. The ANITA observation 
of UHECRs of energies around 10 19 eV at distances beyond 100 km provides experimental support 
for this hypothesis [8]. Next we consider that the field strength scales linearly with energy [9], and 
that the lateral profile drops exponentially e(a) = KEexp(— Ltana/A), where a is the angle under 
which the shower maximum is seen from the antenna position with respect to the shower axis, K 
a linear coefficient, A the lateral profile exponential attenuation parameter and L the projection on 
the shower axis of the distance between the antenna and shower maximum. We define a detection 
threshold angle (%, below which the field is strong enough for detection. Experimental results 2 
indicate a value a t h ~ 3° fori? = 10 17 eV. We can then write tan(a t h) = £ log(f?/10 17 ) + 0.05. 

We consider that a shower is detected if its energy exceeds a value f? t h and if a cluster of 8 
antennas or more are in direct view of the t decay point, at a distance between 20 and 120 km, 
and inside a cone of half angle a t h around the shower axis. We evaluate two different approaches: 
conservative (E t h = 10 17 eV, A = 200 m, yielding a t h = 7° at 10 18 eV and 12° at 10 19 eV) and 
aggressive (f? t h = 3 x 10 16 eV, A = 400 m, yielding a t h = 12° at 10 18 eV and 20° at 10 19 eV). 

Note that the Cherenkov effect is not taken into account in this study. This will be done in the 
up-coming full Monte-Carlo study of the GRAND neutrino sensitivity. 

Results and implications. Assuming a 3-year observation on an area of 60000 km 2 with no 
neutrino candidate, a 90% C.L. integral limit of 6.6 (3.1) x 10 10 GcV 1 cm~ 2 s 1 can be derived 
for an E~ 2 neutrino flux in our conservative (aggressive) scenario. This is a factor > 5 better 
sensitivity than that of other projected giant neutrino telescopes for EeV energies [11]. 

This preliminary analysis also demonstrates that mountains constitute a sizable target for neu¬ 
trinos, with ~50% of down-going events, corresponding to neutrinos interacting inside the moun¬ 
tains (Fig. 1). It also appears that specific parts of the array (large mountains slopes facing another 
mountain range at distances of 30 — 80 km) are associated with a detection rate well above the aver¬ 
age. We thus realized that a factor of 10 improvement in sensitivity - corresponding to cosmogenic 
neutrinos event rates between 10 to 100 per year - could be reached with a factor of 3 increase 
in the detector area, provided the detector is composed of several sub-arrays of smaller size (few 
10000 km 2 ) deployed solely on favorable sites. This is the envisioned GRAND setup. 

The expected angular resolution on the arrival direction of the events detected was computed 

^Vertical showers of energy ~ 10 17 eV are detected at distances up to ~300 m from shower axis [9, 10]. 
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Figure 1: Left: Effective area of the 60 000 km 2 simulated set-up as a function of zenith angle 6 for various 
initial v T energies. Zenith angles < 90° correspond to down-going trajectories. Right: Flight distance 
probability density functions for ts produced by the interaction of v T of various initial energies in a 30 km 
thick wall. For both plots black: 10 17 5 eV, red: 10 18 5 eV, green: 10 19 5 eV, blue: 10 205 eV). 


analytically following [12]. A mean resolution of 0.05° should be achievable for a 3 ns precision 
on the antenna trigger timing. Energy estimation will be challenging, mainly because the r decay 
point is unknown. However, the v T energy is correlated to measurable parameters (Fig. 1) which 
could provide a handle to statistically create the neutrino energy spectrum. 

3. Background rejection 

A few to a hundred cosmic neutrinos per year are expected in GRAND. The rejection of events 
initiated by high energy particles other than cosmic neutrinos should be manageable, as: i) the flux 
of atmospheric neutrinos is negligible above 10 16 eV [13], ii) the rate of showers generated by muon 
decay above GRAND is expected to be few per century [14], and iii) rejection of all trajectories 
reconstructed down to 1° below the horizon would strongly suppress the measured rate of EAS 
initiated by standard cosmic rays, without significantly affecting GRAND neutrino sensitivity. 

The event rates associated to terrestrial sources (human activities, thunderstorms, ...) are diffi¬ 
cult to evaluate, but a conservative estimate can be derived from the results of the Tianshan Radio 
Experiment for Neutrino Detection (TREND). TREND [12] is an array of 50 self-triggered anten¬ 
nas deployed in a populated valley of the Tianshan mountains, with antenna design and sensitivity 
similar to what is foreseen for GRAND. The observed rate of events triggering antennas over a sur¬ 
face > 1 km 2 was 15 events/day in TREND, which scales to a safe estimate of ~ 10 9 events/year for 
a 200 000 km 2 array. A background rejection rate better than 10 9 is therefore probably necessary to 
reach the expected sensitivity of GRAND. Meeting this requirement without affecting the neutrino 
detection efficiency is a key challenge for the GRAND project. 

Amplitude patterns on the ground (emission beamed along the shower axis or signal enhance¬ 
ment along the Cherenkov ring 3 ), as well as wave polarization [15] are strong signatures of v- 
initiated air showers that could provide efficient discrimination tools. These options are being in¬ 
vestigated within GRAND, through simulations and experimental work. In 2016 the GRANDproto 

3 https ://indico.in2p3 .fr/event/10976/session/3/contribution/40/material/slides/0.pdf 
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project [16] in particular, will deploy an hybrid detector composed of 35 3-arm antennas (allowing 
a complete measurement of the wave polarization) and 24 scintillators, that will cross-check the 
EAS nature of radio-events selected from a polarization signature compatible with EAS. 

A project is also being initiated to perform similar tests at Auger-AERA [17]. Its goal is to 
make use of the good detection efficiency of the Auger fluorescence detectors down to large zenith 
angles to complement the GRANDproto study with very inclined EAS trajectories. 

4. GRAND engineering array 

Detailed design of the GRAND detector will be refined once the technical requests associated 
to the definite science case of GRAND have been defined. It is however certain that the GRAND 
detector will be a technical challenge given its unprecedented size. Its huge number of units neces¬ 
sarily implies that they should be fully autonomous over periods of years, with a very low failure 
rate. The driving principles of GRAND to ensure its success will be to stay as modular and basic 
as possible and rely on standard and validated industrial-scale solutions whenever possible. Be¬ 
fore considering the complete GRAND layout, an engineering array of size 1000 km 2 will be 
deployed in order to test the proposed technological solutions. This array will obviously be too 
small to perform a neutrino search, but cosmic rays should be detected above 10 18 eV. Their recon¬ 
structed properties (energy spectrum, directions of arrival, nature of the primary) will enable us to 
validate this stage, if found to be compatible with the expectations. The absence of events below 
the horizon will also test our EAS identification strategy. 

5. High-energy neutrino astronomy with GRAND 

UHECRs are likely produced in extragalactic sources, given the strength of Galactic magnetic 
fields and the lack of correlations with the Galactic plane. Some fraction of their energy is con¬ 
verted to high-energy neutrinos through the decay of charged pions produced by interactions with 
ambient matter and radiation. This can happen in the source environment or during the flight in the 
intergalactic medium (cosmogenic neutrinos). The range of expected cosmogenic neutrino fluxes 
can be calculated precisely, and depends mostly on parameters inherent to the cosmic rays (Fig. 2) 
[18]. The calculated level of high-energy neutrino fluxes produced at the sources depends on the 
modeling of the acceleration region. But robust flux estimates have recently been derived for some 
sources [19, 20, 21, 22], arguing that many scenarios can be tested with increased sensitivities. 

The sensitivity of GRAND should guarantee the detection of cosmogenic EeV neutrinos. For 
reasonable source scenarios (Fig. 2), GRAND aims at collecting in the order of 100 events per year 
above 3 x 10 16 eV, which helps to understand the underlying components. An energy resolution of 
50% would enable us to picture the shape of the diffuse energy spectrum, but the most important 
information would stem from the precise (< 0.1°) angular resolution on the arrival directions. 
Cross-correlating the position of sources and of events, it will be possible to discriminate between 
cosmogenic neutrinos and those produced in the source environment. The former should indeed 
correlate with the large-scale structures at high redshift where the integrated flux is maximal. 

The ideal way to identify high-energy neutrino sources would be to observe a point source 
—a transient source is a likely candidate due to the power and density constraints from UHECR 
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Figure 2: Left: Differential sensitivity limit of the 60000 km 2 simulated setup (red dashed lines, thin: 
conservative, thick: aggressive), of the projected 3 times larger GRAND array (red thick long-dashed), and 
for other instruments, and estimated theoretical cosmogenic neutrino fluxes for all flavors [18], The blue line 
gives the most pessimistic fluxes (pure iron UHECR injection and low maximum acceleration), while the 
gray shaded-region indicates the “reasonable” parameter range. Right: Integrated exposure of the simulated 
setup over 3 years at 10 18 5 eV, in units of cm 2 s, in Galactic coordinates. 


observations (e.g., [23]). GRAND opens this possibility with its excellent spatial resolution and 
its sky coverage. In the simulation layout, the instantaneous field of view of GRAND is a band 
corresponding to zenith angles 80° < 9 < 100°, and its integrated exposure 10 16 cm 2 s for 3 
years in the energy range ~ I0 17 20 cV, over a large portion of the sky (Fig. 2). Short-lived (< 
day) transients have a low probability of being spotted, but for longer transients (blazar flares, 
ultraluminous supernovae, magnetars etc.), post-analysis at the location of existing transients and 
stacking searches can be done efficiently. Depending on the background discrimination efficiency, 
GRAND would be able to send alerts to other experiments via AMON-type structures [24], 

The design of GRAND is still preliminary, and can be adapted to accommodate further re¬ 
quests related to the high-energy neutrino astronomy. The large array size and the frequency range 
of GRAND - if extended to 300 MHz - could be profitable for other physics goals (such as the 
study of the Epoch of Reionization [25] or the search for Fast Radio Bursts [26]). We are presently 
assessing the technical cost of pushing the detector towards these directions. 

6. Ultrahigh Energy Cosmic Rays 

The energy spectrum of cosmic rays shows an increased steepening at around 10 19 6 eV [27, 
28]. The origin of this steepening is still unclear. The Telescope Array (TA) analysis is consistent 
with a dominance of protons up to the highest energies [28], while the Auger experiment favours 
a mixed composition [29] . Also TA has reported an anisotropy in the arrival directions of UHECR 
in the Northern hemisphere [30], whereas the Auger collaboration has not identified a significant 
anisotropy at the highest energies in the South [31]. Due to the absence of photons at high energies, 
top-down and other exotic scenarios of UHECR production are disfavoured. 

Given its effective area, GRAND would be the largest UHECR observatory on ground, and 
should provide valuable information to help resolve these puzzles if precision of reconstructed 
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parameters is satisfactory. Determination of the UHECR composition in particular is a key issue. 
Detailed simulation studies will be performed to estimate performances reachable by GRAND. 

7. Conclusion 

The GRAND project aims at building a next-generation neutrino telescope composed of radio 
antennas deployed over ~ 200000 km 2 . Preliminary simulations indicate that a sensitivity guaran¬ 
teeing the detection of cosmogenic neutrinos could be achievable. GRAND would also be a pow¬ 
erful instrument for UHECR detection at energy > 10 19 6 eV with high statistics, and possibly for 
other science objectives (Epoch of Reionzation, Fast Radio Bursts). Work is ongoing to determine 
precisely the achievable scientific goals and the corresponding technical constraints. Background 
rejection strategies and technological options are also being investigated. 
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